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Abstract. The quasispecies theory is studied for dynamic replication 
landscapes. A meaningful asymptotic quasispecies is defined for periodic 
time dependencies. The quasispecies' composition is constantly changing 
over the oscillation period. The error threshold moves towards the posi- 
tion of the time averaged landscape for high oscillation frequencies and 
follows the landscape closely for low oscillation frequencies. 

The quasispecies theory, put forward by Eigen in 1971 H, and subsequently 
studied by Eigen, Schuster, McCaskill and coworkers ||, HH, i s nowadays one 
of the classical theories of self-replicating entities. Its prediction of an error 
threshold, above which the self-replication ceases to produce useful offspring, 
has important implications for the origin of life. The error threshold effectively 
limits the amount of information the entities can carry, thus placing an up- 
per bound on the complexity self-reproducing information carriers can achieve 
without sophisticated error correction mechanisms. 

Although completely static environments are unrealistic in any case apart 
from experiments in perfectly controlled flow reactors, the quasispecies theory 
has so far been considered mainly in static replication landscapes. Neverthe- 
less, even under fixed environmental conditions can the replication rates of RNA 
molecules, for example, change because of changing concentrations of template 
and replica f| . Jones |(], |t]] has studied underlying time-dependencies which are 
identical for all sequences. Contrasting to that, we want to focus on replication 
landscapes with individual time-dependency for each sequence. One of the rea- 
sons for the neglect of individually changing replication coefficients in earlier 
work is probably the fact that for arbitrary temporal changes an asymptotic 
quasispecies cannot be defined. However, a meaningful definition is at hand for 
time-periodic replication landscapes, as we are going to show below. 

We start from the discretized form of Eigen's evolution equation B, lin- 
earized with the appropriate transformation B [lO| . Due to space limitations, we 
cannot repeat the arguments leading to that equation here. For details about 
this calculation, the reader is referred to j|. We use the same notations as are 
used there. Additionally, we define the error rate R = 1 — q, which gives the 
probability that a single symbol is copied erroneously. The string length will be 
denoted by I throughout this paper. 



The vector of the unnormalized sequence concentrations y(t) evolves accord- 
ing to 



y(t + At) = AtW(t) + 1 y(t). 



(1) 



Here, W(t) is the replication matrix W(i) = QA(i) — D(t). We assume the matrix 
W(i) is periodic with period T = nAt, n £ N, with some chosen discretization 
time step At < T. After iteration of Eq. @), we obtain for t' = t + (At, (eN 



y(t') 




AtW(t + vAt) + l } y{t). 



(2) 



where T{.} stresses that the product has to be evaluated in the time order given 
by the iteration. With the definition of the matrix (or operator) 



X :=T 



Yl [AtW(uAt) + l 



(3) 



which maps y(0) onto y(T), we are now able to write down the solution of the 
discretized differential equation Eq. (|]) for the initial condition y(0) as 



vit) 




X m y(0) . 



(4) 



where the time t has been subdivided into t — mT + (At, with £ < n and 
m, C G N. 

If we observe the system in time steps of the period length T, the system 
appears to evolve in a static replication landscape, which is defined by X. The 
asymptotic steady state for the oscillation phase £ = is therefore given by 
the normalized Perron eigenvector <p Q of X For < ( < n, the steady 
states are found by application of T{. . . } from Eq. (Q) to (fi and subsequent 
normalization. 

Let us now study quantitatively the effects a periodic replication landscape 
has on the prominent quasispecies. As the first step into that direction, we start 
from the Swetina-Schuster landscape and introduce small oscillations in the 
master sequence's replication coefficient Aq. For reasons of simplicity, we set all 
decay constants equal Di(t) = D, because then they drop out of Eq. (|l|) during 
the foregoing linearization. 

We will write the time-dependent replication coefficient Ao(t) as 



A (t) = A ,sexp[e/(i)] 



(5) 



where Aq^ is the replication coefficient in the static landscape, f(t) is a T- 
periodic function and e is the oscillation amplitude. For e = the corresponding- 
static landscape is reached. The other replication coefficients are equal A% = 
■ ■ ■ = A\ = A and constant. We will choose A so small that the condition A <C 
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Fig. 1. The steady state oscillations of the master sequence in cyclically changing 
environments with different oscillation periods T. Parameters used are I = 2, 
A 0<s = e 2 - 4 , A 1 (t)=A 2 (t) = 1, e = 0.2. 

Ao(t) is satisfied for all t and e«l. This assures that we see a clear transition 
from the static case to the dynamic case, and additionally, that the changes in 
the master sequence's abundance can be directly related to the changes in Ao(t). 
One of the simplest forms the function f(t) in Eq. (|J) can take on is 

f(t) = sin(wi) with w = 2ir/T. (6) 

In the following, we will shortly discuss the influences of different frequencies u> 
and amplitudes e for this time dependency. 

As response to the oscillation of the replication coefficients, a modified os- 
cillation is found in the concentration xq of the Master sequence (Fig |l|). For 
increasing frequency uj, the amplitude of the xq oscillation decreases and a phase 
shift strengthens. This behaviour is due to the finite time a reaction system as 
described by Eigen's equation needs to settle into equilibrium. In constant envi- 
ronments, the asymptotic species distribution is approached in exponential time, 
with the relaxation time scale r set by the difference between the largest and the 
second-largest eigenvalue of W. For the oscillating environments the relaxation 
time needs to be compared to the period T. If T ^> r, the system is virtually in 
equilibrium for arbitrary (asymptotic) times t, whereas for T ~ r the changes 
cannot be tracked anymore and phase shift as well as amplitude damping of the 
response sets in. For T«t the response amplitude gets fully damped and the 
system gets identical to one with the time-averaged replication coefficients. Inter- 
preting the Aq and xq time-dependence as input and output signal, the system 
acts as a low pass analog filter, in analogy to observations made in population 
genetics models with dynamic fitness landscapes [[l3[ |l4|, [TJ| ^6). Moreover, for 
small e, the filter works linear, which means that a sinusoidal oscillation is found 
in xq , whereas this linearity is quickly destroyed for increasing eQ. 

1 Details on this analysis can be found in |17J|. 





Fig. 2. The quasispecies distribution as a function of the error rate R. Aq,s = 
e 2A ,A = 1, e = 2, T = 100. Two different oscillation phases are shown, left: 
( = n/2, right: ( = 0. 

We will now focus on the influence of a time-dependency as given in Eq. (||) 
onto the error-threshold. In accordance with the above, we have to distinguish 
between different dynamic regimes. For T -C r, a sharp error-transition occurs 
at R^ v , which denotes the error-threshold of a system with time-averaged Aq = 
^ J Ao(t)dt. Contrasting to that, a moving error-transition can be found for 
T ^> r at approximately R* (t) , which denotes for any given t the error-threshold 
in a constant landscape with A = A (t). R*(t) lies between i?* nax and i?^ in , 
which correspond to Aq = max t Aq (t) and Aq — min f Aq (t) , respectively. In the 
intermediate cases T w r, the numerical simulations (see Fig. ||) show that the 
error-threshold R*(t) oscillates within a smaller interval than [i?^ in , i?j^ ax ]. 

These findings allow to draw a phase diagram as displayed in Fig. [| For 
low T, we observe the standard separation into an ordered phase (below the error 
threshold) and a disordered phase (above the error threshold). With increasing 
oscillation period T, a third, new phase appears between the two. In this phase, 
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Fig. 3. Schematic phase diagram for a time-dependency like Eq. 



we observe — for a fixed error rate R — an alternation between a fully developed 
quasispecies and a completely disordered system. The population seems to be 
moving back and forth over the error threshold. Therefore, we call this new phase 
the temporarily ordered phase. Since a similar phase diagram can be expected 
for any periodic landscape with finite i?* v and i?* lin ^ -Rm ax , we believe that 
such observations could also be made in typical AL simulations such as Tierra or 
Avida ]lS[ ] provided with the appropriate replication landscape. The temporarily 
ordered phase would for a finite population in a rugged landscape have the effect 
of causing a random drift over the landscape at some times and a localization 
around a local master sequence at other times. 

Upon completing this work, we became aware of Ref |]l9| , in which a different 
approach towards dynamic replication landscapes is given, using a stochastic 
time dependency in the landscape. The results presented there cannot directly 
be related to our findings here, because the equivalent to R* av vanishes in fl9|j , 
while the equivalents to R^a* and -R^in take on the same finite value. This leads 
to a different phase diagram than the one we observe here. 
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